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Abstract. Some selected topics in photoproduction are reviewed. In particular, the 
focus is on quasi-real photon-proton scattering as measured at HERA. Single & dijet 
rates and jet profiles are discussed, as are open charm production and inelastic char- 
monium production. The many intriguing theoretical issues which can be addressed 
through study of these processes are highlighted. 

INTRODUCTION 

In this talk, I’m going to present a snmmary of recent developments in pho¬ 
toproduction at HERA. By “photoproduction” I mean yp interactions in which 
the photon is close to being on its mass shell. Furthermore, I shall focus on pro¬ 
cesses in which there exists at least one hard scale - so that we have the chance 
to use perturbation theory. I’m going to divide the talk into 4 parts: jets (single 
& dijet rates); jet prohles; open charm production (Zi)*’s); inelastic charmonium 
production. I won’t talk about processes with rapidity gaps, since these diffractive 
processes are covered in other people’s presentations. 

I JET CROSS-SECTIONS 

I wish to start by recalling some basic terminology. It is convenient to talk 
about two classes of photoproduction event: those in which all of the energy of 
the incoming photon is delivered to the hard subprocess (these are called “direct” 
processes) and those in which only part of the incoming photon energy goes into 
producing the hard subprocess (these are called “resolved” processes). 

Fig. 1(a) shows the lowest order contribution from these two types of process in the 
case of jet production. The dotted lines show where we chose to separate the hard 
process from the non-perturbative parton distribution functions. The factorization 
scales and Pp are in principle arbitrary but a good choice would be to pick them 
to be ~ pt (i.e. the jet transverse momentum) since this choice will account for 
the fact that the incoming partons have a large region of transverse phase space 
into which they can radiate other partons. Nonetheless, it is certainly true that 
the theoretical calculation will exhibit a strong sensitivity to any variation in the 
factorization scales since the 2 —2 hard subprocesses contain no dependence upon 
the factorization scales whereas the parton distribution functions do. For example. 
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FIGURE 1. LO (a) and NLO (b) direct and resolved contributions to jet photoproduction 

we could effectively turn off all parton evolution by picking the factorization scales 
to be much smaller than pt- At this level of calculation, it is clear what we mean 
by “direct” and “resolved” processes. 

In Fig.1(b), the 0{as) corrections to the jet production process are shown. We 
have to account for the fact that the hard subprocess can radiate off 3 partons into 
the hnal state as well as the virtual corrections to the 2 parton hnal state. Again, 
we can classify diagrams as being direct or resolved. However, there are collinear 
divergences associated with the 3 parton hnal state. This means that the 3 parton 
hard subprocesses are dependent upon the factorization scales. Let us focus on the 
dependence. By varying we vary the amount of the next-to-leading order 
(NLO) 2 —3 processes that should really be attributed to the lowest order (LO) 
resolved process. In particular, the amount of NLO direct that should rightfully be 
factorized into the LO resolved contribution depends upon the choice of In this 
way, the strong dependence upon the factorization scales is reduced by including 
the NLO contribution. Also, we can no longer uniquely dehne the separation of 
direct and resolved processes. 

A Single jet inclusive 

NLO calculations of the single jet inclusive rate have been available for a while 
now [1]. Comparison with the data [2-4] should tell us about the parton content 
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FIGURE 2. Inclusive single jet cross-section integrated above Et = 8 GeV (a), 11 GeV (b) and 
17 GeV (c). The data [4] are compared with the NLO calculation of [5]. Dotted line shows the 
direct contribution, dashed line the resolved and solid the total. 


of the photon, and in particular its gluonic content (which is not really constrained 
from elsewhere). However, in just the region where one expects sensitivity to the 
gluon in the photon the NLO theory is not working! Forward jets (i.e. those heading 
largely in the proton direction and hence having positive pseudo-rapidity, if) with 
modest values of transverse energy, Et, are most sensitive to the gluon content of 
the photon. Fig.2 compares the 1993 ZEUS data with the theoretical calculations 
of Klasen, Kramer & Salesch [5] which use the GS(HO) photon parton densities 
[6] and the MRSD- proton densities [7]. Comparison reveals that the theory is 
well below the data for forward jets at the lower Et values (this conclusion is also 
supported by the HI data and the ZEUS 1994 data). The problem seems not so 
bad as Et is increased. One might argue that the difference between theory and 
data should be attributed to the fact that the photon parton density functions need 
adjusting, i.e. more gluons are required at low x. However, one should be cautious. 
These forward jets are also much broader than expected from Monte Carlo studies. 

It is possible that the broadening of the forward jets comes about in the presence 
of a large underlying event. Models which attempt to account for some or all of this 
physics [8,9] have received some support from the data [3]. In particular, simple 
models of multiple parton scattering produce large effects in the forward region. 








e.g. see [9]. Multiple parton scattering is anticipated on the grounds that forward 
jets at low Et are produced as a result of interactions between slow partons in 
the colliding particles. We know that QCD predicts a proliferation of these slow 
partons, and as such it may well be that more than one pair of them can interact 
in each yp interaction. 

In summary, the HERA experiments are providing precise data which is ready 
for comparison with the theory. However, before we can safely extract the photon 
parton distribution functions we really need to understand what is going on with 
the forward jets. I think the issue of the forward jets is potentially very interesting. 
If multiple interactions really are present, then we are looking at a new regime of 
QCD. 


B Dijets 


Data on two or more jets [10-12] provides us with further options to test QCD 
and understand the nature of the “strongly interacting” photon. In particular, the 
two highest jets can be used to compute the variable 
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At lowest order, = 1 for direct events and x^ <1 for resolved events. A cut on x^ 
therefore provides a physical dehnition of direct and resolved event classes. ZEUS 
has dehned direct enriched and resolved enriched samples by separating events 
according to a cut at x^ = 0.75. The direct enriched sample is very sensitive to the 
small-x gluon content of the proton: the more backward the dijets, the lower the x 
values in the proton that are probed. Conversely, the resolved enriched sample is 
sensitive to the gluon content in the photon. In addition, NLO calculations for the 
dijet rates are now available for comparison with the data [13,14]. Let’s summarize 
the situation as it stands right now. 

Eor x^ > 0.75, as can be seen from Fig.3, the NLO theory does a good job. 
However (and this is not shown in the plot) there remains qnite a large contamina¬ 
tion from the large-x part of the photon qnark distribntion fnnctions. This arises 
because of the harder form of the photon qnark densities. To unravel the effects 
of the low-x glnons in the proton from the large-x qnarks in the photon reqnires a 
tighter cnt on x^. This is feasible, for example for Et > 30 GeV and 250pb ^ of ep 
data one expects around 4500 events with x^ > 0.9 [17]. Also, comparison between 
data and theory reqnires that the same jet algorithm is used. To facilitate a clean 
comparison between data and theory, the ZEUS collaboration has started to used 
the /crcluster algorithm [18]. 

For x^ < 0.75 Fig.3 reveals that the theory falls well below the data for the lowest 
Et forward dijets. The effect exhibits a strong dependence npon the Et cnt, which 
snggests that it cannot be explained by modifying the parton distribntion fnnctions 
of the photon in any sensible way. Presnmably, this is the same problem as that 




FIGURE 3. Inclusive dijet cross-section for direct and resolved enriched processes as a function 
of the average jet rapidity, fj, and integrated over the rapidity difference, Arj, of the two highest 
Et jets, i.e. |A? 7 | < 0.5, and Et as shown. The LO and NLO theory of [13] is compared to the 
data [10]. The GRV [15] and CTEQ3M [16] NLO MS parton densities were used. 


which we encounter for the single jets. Again, multiple interactions help fix the 
problem. 

The dijet measurements shown in Fig.3 have been made with a cut on the mini¬ 
mum Et of the jets and the cut is the same for both jets. This introduces a further 
theoretical problem. This arises because most of the jets will be produced around 
the minimum allowable Et, i.e. the typical difference between the jet transverse 
momenta, ApT, will be small. So, the 3 parton final state (which is present in the 
NLO calculation) must have one of the partons either collinear with another, or 
very soft. The collinear configuration is easy to deal with (it is factorized) but the 
soft parton emission leads to a InApT contribution. This large logarithm signals 
that multiple soft parton emission is important. These soft parton effects can be 
studied by looking explicitly at the ApT distribution of the dijets or they can be 
avoided by making a cut which keeps away from ApT ~ 0. 

Complementary to the jet rapidity and Et distributions is the distribution in 
the polar angle of the jets, as defined in the dijet centre-of-mass frame. Such a 
measurement is insensitive to the parton density functions but sensitive to the 
subprocess which drives the jet production. Fig.4 shows the data [11] and the good 


















ZEUS 1994 



FIGURE 4. Dijet angular distribution 


agreement with NLO theory [14], The steepening of the cos0* distribution as 
I cos 0*1 —1 is greater in the resolved enriched sample than in the direct enriched 
sample due simply to the fact that the resolved process proceeds predominantly 
through f-channel gluon exchange whilst the direct processes go via t-channel quark 
exchange. 

In conclusion, the dijets provide information which complements the single jet 
measurements. The data is now reaching a high level of precision, and comparison 
with NLO theory has revealed a number of pressing issues. In particular, we need 
to understand better the forward jets, ensure that we are using the most convenient 
jet algorithm, and be aware of any sensitivity from soft parton emission. Once these 
issues have been addressed, we can expect to gain further insight into the gluon 
content of both the photon and proton. 


II JET SHAPES 

Now I want to address the question; “What is a jet?”. An observable which has 
been used to investigate this question is the shape variable ip{r, R, ET,'f]) which is 
the fraction of the jet (of cone size R) energy lying within a cone of size r [19]. The 
lowest non-trivial order for this quantity involves the 3 parton matrix elements. To 
compute to order requires the 4 parton hnal state matrix elements [20]. Klasen 
& Kramer have been able to £t the ZEUS data on jet shapes [21] but only after 
introducing the parameter, Rsep, into their calculations [22]. Let’s take a closer 
look at the i?sep parameter. 

It’s a good idea to first remind ourselves of the steps involved in building jets using 
a cone algorithm [23]. The ZEUS jets are found using an iterative cone algorithm 
(PUCELL). This algorithm defines some threshold energy, Eq, and calorimeter cells 
above this energy define the seeds for the jet finding. A jet is formed by summing 





all cells within a distance i? (in 77 — 0 space) of the seed cell. If, after this, the jet 
axis no longer coincides with that of the original seed cell then the new jet axis is 
used and all cells within R of it are used to redehne the jet. This process is iterated 
until a stable jet direction is achieved. It is possible that some of the jets overlap 
and so some criterion needs to be established to decide what should be done with 
these jets. For example, if one of the overlapping jets has more than 75% of its 
energy in common with another jet then one could merge the two jets, whilst if the 
overlap energy is less then one could split the two jets by assigning particles to the 
jet to which they are nearest. Finally, jets are accepted if their Et is larger than 
some minimum value. In the theoretical case, two partons are combined if they he 
within a distance Rij of each other, where 


Rij < Klin 


{Exi + Exj)R 
max{ET^,ET,y • 


( 2 ) 


We have introduced i?sep into this dehnition, it is introduced so as to simulate the 
effect of seed hnding. For example, two partons with equal ERs lying a distance 
2R apart could be considered to form a single jet (the partons lying on the edge of 
the jet cone). Such a conhguration could not occur in the experimental algorithm, 
since there is no seed cell in between the two partons. As a result, for a 3 (or less) 
parton hnal state, picking i?sep = R should match the theoretical and experimental 
dehnition of the jets. Notice that the 3 parton hnal state does not call for any jet 
merging. By htting i?sep to the data, one is accounting for higher order ehects by 
htting a single parameter (to ht the data Klasen & Kramer need to vary Rsep from 
1.3 to 1.8 for jets with 17 ^ —1 to 17 ~ 2 respectively). For example, parton shower 
and hadronization ehects (plus, maybe, multiple parton scattering) all ehect the 
hnal state and we do not really learn much about them by tuning the theoretical 
calculation. Fig.5 demonstrates that forward jets are broader, as I claimed earlier. 
It also shows how multiple parton interactions can help improve the situation. It 
really is a challenge to theorists to explain the variation of Rsep- 

Steps in this direction require a truly NLO calculation of the jet shape variable. 
As mentioned above, this means we need to consider the 4 parton hnal state. At 
this point, we realise that the iterative cone algorithm is not infra-red safe. By 
this we mean that it is not insensitive to the addition of a soft parton [20]. To see 
this consider the 4 parton hnal state illustrated in Fig. 6 . The middle of the upper 
3 arrows denotes a soft parton, which has energy just above the threshold energy, 
Ep. The presence of this seed in between the two hard partons would cause all 3 
partons to be merged into a single jet. So our shape variable is sensitive to the 
cell threshold energy. This is an acute problem at any hxed order in perturbation 
theory (beyond LO). However, a strong Eq dependence isn’t seen in Monte Carlo 
studies. Let me outline the basic idea as to why this should be so [20,24]. Denote 
the probability for emitting a single parton into the overlap region. Pi, where 
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FIGURE 5. Jet shape from single jet production. Comparison between data [21] and PYTHIA, 
both with and without multiple parton scattering (MI), is shown. 
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FIGURE 6. A four parton final state 














This probability exponentiates on summing over all possible numbers of soft par- 
tons, i.e. 

Paii~l-e-^b (4) 

So as Po —^ 0 a soft parton is guaranteed to be in the overlap region, i.e. Pali —^ 1- 
Thus we see that, on expanding the exponential, any hnite order in is sensitive 
to large IuPq terms but that this sensitivity goes away on summing to all orders 
(as is done in the parton shower routines). As pointed out by Steve Ellis [25] this 
problem can be fixed simply be redehning the jet algorithm, so that the mid-point 
between jets is always treated as a seed, regardless of any activity. Alternatively, 
the kt cluster algorithm avoids such problems. I think it is fair to say that studies 
on jet structure are reaching the point where theory really is starting to address 
the important issues. We should learn much over the coming few years. 

Ill OPEN CHARM 

The photoproduction of charm quarks at large pt is a process which involves two 
large scales, pT,'mc ^ A and as such, makes life more complicated from the theo¬ 
retical point of view. Good data, which can be expected in the future (especially 
if the charm can be tagged using a microvertex detector), will surely shed light 
on this intriguing area. At present, there are two main routes used in theoretical 
calculations. I’ll start by describing briefly each. 

Massive Charm: The charm quark mass is considered to provide the hard scale, 
as such charm only ever appears in the hard subprocess and there is no notion of 
radiatively generated charm in the sense of parton evolution [26]. This means that 
terms ~ \n{p‘^/ml) are not summed to all orders (in the parton distribution and 

fragmentation functions). As such, we might expect this approach to become less 
accurate when pf. 3> However, it does provide a systematic way of accounting 
for charm quark mass effects, which will be important for rric ~ Pt- 

Massless charm: In this approach, the charm quark is treated as massless 
(above threshold), and as such is treated like any other light quark in the parton 
evolution equations and hard subprocess cross-sections. The In(p^/mc) terms are 
now summed to all orders, but charm quark mass effects are ignored. So, this 
approach should get better as pr/^c increases. 

At present the data [27-29] lies in the intermediate region where pt < 10 GeV, 
i.e. it is not clear which, if any, of the two approaches should be used. In order to 
compute the inclusive D* rate, one needs the appropriate fragmentation function. 
Either the Peterson [30] form or the x°‘{\ — form do a good job, and have been 
well constrained by data. 

Initial comparison between theory and data suggested that the massive charm 
calculation was too low, e.g see [29]. More recently, the situation has changed some¬ 
what [31,32]. The full NLO calculations require that the fragmentation function 
be consistently extracted from the data. When this is done, Gacciari and 
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FIGURE 7. Comparison of the D* cross-section data [27] with the calculation of [31]. 


Greco find that the theoretical predictions are increased signihcantly relative to 
what is found using the softer (LO) fragmentation functions [31]. This is true for 
both massless and massive charm calculations and, within theoretical uncertainties, 
both are now consistent with the present data. For example. Fig.7 is a plot from 
[31], which shows that the NLO theory assuming massless charm agrees well with 
the data. 

More data at large pt and increased statistics at intermediate pt will certainly 
help in our study of the interplay of \n{p‘^/ml) and ml/p^ effects. In addition, for 
Pt 3> rric we have the possibility to study the “intrinsic” charm within the photon 
(charm in dijets offers good prospects here). 


IV CHARMONIUM 

Originally, inelastic photoproduction of charmonium, e.g. J/'ip, was advertised 
as an ideal way to extract the gluon density in the proton (since it is driven by 
photon-gluon fusion into a QQ pair). More recently, NLO calculations have put a 
dampener on this goal (as we shall see shortly). However, there has been a great 
deal of recent interest in the non-relativistic QCD (NRQCD) approach to heavy 
quarkonium production, and the inelastic photoproduction of heavy quarkonia pro¬ 
vides the ideal opportunity to test NRQCD. 

Bodwin, Braaten and Lepage derived a factorization formula which describes 
the inclusive production (and decay) of a heavy quarkonium state [33]. In the 
case of photoproduction. Fig.8 shows the lowest order contribution. The NRQCD 
factorization formula for the corresponding cross-section reads 




H 



FIGURE 8. Leading contribution to quarkonium photoproduction. 

da{H + X) = Y,da{QQ[n]+X){0^). (5) 


X denotes that the process is inclusive, dd{QQ[n] + X) is the perturbatively cal¬ 
culable cross-section for yp — QQ + X and it can be written as a series expansion 
in aslmg). The QQ pair is produced with quantum numbers n. The matrix ele¬ 
ment, (O^), contains the long distance physics associated with the formation of the 
quarkonium state H from the QQ state - it is essentially the probability that the 
pointlike QQ pair forms H inclusively. The typical scale associated with this part 
of the process is ~ mqv which is much smaller than mg {v is the relative velocity 
of the QQ pair, and is small for heavy enough quarks). This hierarchy of scales 
underlies the NRQCD factorization. Note that the QQ state is not restricted to 
having the same quantum numbers as the meson. Fortunately, there exist “velocity 
scaling rules” which allow us to identify which states, n, are the most important. 
More precisely, the “velocity scaling rules” order the operators (O^) according to 
how many powers of v they contain, i.e. relativistic corrections can be computed 
systematically. 

The NRQCD approach therefore provides us with a systematic way of computing 
inclusive heavy quarkonium production (modulo corrections which are suppressed 
by powers of ~ A/mg). The strategy is hrst to organise the sum over n into an 
expansion in v and then to systematically compute dd order-by-order in asirng). 
Technically, we do not a priori know where our efforts are best placed, i.e. do we 
work at lowest order in v and to NLO in as or do we attempt to work at higher 
orders in v, but computing each hard subprocess to lowest order? We need to know 
V in order to judge better what to do. 

One dual word before moving on to discuss J/'ip photoproduction. For small pt, 
NRQCD factorization is likely to break down, due to contamination from higher 
twist effects. Also, one expects breakdown of the NRQCD approach in the vicinity 
of the elastic scattering region, i.e. z ^ 1 where 2 : is the fraction of the photon 
energy carried by the quarkonium (see later). 

Inelastic photoproduction of J/ip is something which has already been measured 
at HERA. Let’s see how the theory shapes up. To lowest order in the velocity 
expansion, [n] = [l,^S'i]. The hrst entry in the square brackets tells us that the cc 



is in a colour singlet state, whilst the second entry tells us the spin and angular 
momentum of the state. Not surprisingly, to lowest order in the velocity expansion, 
the cc must be produced with the same quantum numbers as the J/'^. This is 
just the colour singlet model (CSM) of old. The lowest order diagram which can 
contribute is shown in Fig. 8 and 


where 0(0) is the wavefunction at the origin (it can be extracted from the electronic 
width of the 7/0). NLO(q;s) corrections have been computed [34] and shown to be 
large. Fig.9 shows that the NLO corrections enhance the LO calculation and lead 
to a reduced sensitivity to the gluon density in the proton. 



FIGURE 9. LO and NLO colour-singlet predictions [35] for the total inelastic J/ip cross-section 
for different parameterizations of the proton parton distribution functions compared to the HERA 
data [36,37] 


One might well ask as to the significance of the resolved photon contribution. It 
is important at small enough [38-40]. In addition, for mesons produced at high 










enough px we have an additional scale to consider and terms which are leading 
in as can be suppressed by powers of ~ This is true for example of the 

diagram shown in Fig .8 relative to that shown in Fig.10. The latter fragmentation 
contribution is higher order in ag, however there is one less hard quark propagator 
and so it will dominate for large enough px- Fragmentation contributions and 
resolved photon contributions are not important in computing the total rate for 
> 0.4 (which is essentially where the data is), and so we expect the curves in 
Fig .8 to be reliable. 



FIGURE 10. Fragmentation contribution to the production of the J/ip- 

Going to NLO in the velocity expansion means moving away from the GSM. For 
the hrst time we encounter colour octet contributions. In particular, the LO(q;s) 
diagram is again that of Fig. 8 ^ now the cc can be formed in one of 5 states, i.e. 
[n] = [8,^S'o], [8,^S'i], [ 8 ,^Po,i, 2 ]- The price one pays for having to convert this state 
into the J/ijj is an extra power of relative to the colour singlet matrix element, 
i.e. 

It turns out that the ~ 0.01 suppression of the long distance matrix elements is 
partially compensated for by a strong enhancement of the corresponding short dis¬ 
tance cross-section. In particular, this is so for the [8,^5'o] and [ 8 ,^Po, 2 ] states. The 
colour singlet matrix element can be extracted from data, e.g. the leptonic width 
of the J/ijj- Likewise, we need to £t these new matrix elements to data (or extract 
them from lattice calculations). It is therefore clear that a test of the NRQD frame¬ 
work requires data from different sources - the challenge being to find a consistent 
description. This is a particularly topical issue, since an explanation of the Teva- 
tron excess of direct J/ijj and tjj' production needs, in addition to fragmentation 
contributions, colour octet contributions [41]. One can use the Tevatron data to fit 
the relevant matrix elements. The validity of this explanation can be confirmed (or 
not) on comparing to data which can be obtained from HERA. Unfortunately, the 
matrix elements which are important at the Tevatron are not so important in J/'ip 
photoproduction for > 0.4. However, the key matrix element for the Tevatron 
([ 8 ,^ S'!]) does play a key role in the region of lower ; 2 , where (large enough px) 


There is a lower order contribution in which no gluon is radiated off, however this would give 
a contribution only at z = 1 which lies outside our region of interest. 



the dominant contribution comes from the fragmentation mechanism via resolved 
photons [39,40]. Another process which is sensitive to the [8,^ S'!] state is the pho¬ 
toproduction of J/tjj + 'y (where the photon is produced in the hard subprocess, i.e. 
not via the radiative decay of a P-wave quarkonium) [42], So, with the anticipated 
increase in statistics, we can really expect to test NRQCD at HERA. Going back to 
the J/'^, there are some weak constraints on the important matrix elements from 
the Tevatron data and these have been used in the theoretical calculations of [43]. 
In Fig. 11 the HERA data on the distribution are shown and seen to compare very 
well with the colour singlet calculation. The colour octet contribution however, is 
much too large at large ; 2 . The shaded band on the theory prediction denotes the 
type of uncertainty expected from the £t to the matrix elements, although this error 
is hard to ascertain with conhdence. Thus the HERA data is not supporting a large 
colour octet contribution at large ; 2 . However, one must be careful in interpreting 
this as evidence against the NRQCD approach, since the z ^ 1 region is sensitive 
to higher order non-perturbative contributions which lead to the breakdown of the 
NRQCD expansion [44]. 



FIGURE 11. Colour singlet and colour octet contributions to the J/b energy distribution [35] 
compared with the HERA data [36,37]. 


In addition to the processes just discussed, increased statistics will allow mea¬ 
surement of other meson states, e.g. %Ij' and T which will certainly further test our 
understanding of QCD. 

Before hnishing. I’d like to mention an alternative to the NRQCD approach which 
is the colour evaporation model, e.g. see [45] (essentially a resurrection of the old 
duality calculations [46]). The production rate is calculated according to 


1 /•2mD daicc) 

a = - dm —;- 

9 J 2 mc dm 


( 6 ) 


where m is the invariant mass of the produced cc pair. This is the cross-section for 





production of charmonium, i.e. it is the sum of production rates for J/'^, 'ip', x ^-nd 
rjc- To get the rate for any particular state, one must multiply by some (empirically 
determined) factor. The factor of 1/9 comes from arguing that the colour octet 
cc convert with probability 1/9 into a singlet state (which must be charmonium in 
this mass range). The colour bleaching is supposed to occur as a result of some 
soft dual state interactions, the precise nature of which is not important. For more 
details on the phenomenological success of the colour evaporation approach, I refer 
to [45] and references therein. 


V FOR THE FUTURE... 

The future bodes well, not least as a consequence of the anticipated increase in 
statistics. Higher statistics will mean more precise studies of open charm production 
over a wide kinematic range (which is vital if we are to develop our understanding) 
and will also allow us to really test our understanding of charmonium production 
and, in particular, NRQCD. The intriguing issue of forward jet production needs 
further work before we can claim to really understand what is happening. NLO 
theory facilitates more precise tests of single and dijet rates using what is becoming 
very good data - it is true to say that we are learning how to really test the theory. 
I didn’t mention the structure of the virtual photon or prompt photon production. 
Good data is starting to accumulate on both these processes and we can anticipate 
that much will be learnt over the next few years. 
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